Mon. Not. R. Astron. Soc. OOP. [TlfT6l C2008') Printed 26 Juno 2008 (MN WT^ style file v2.2) 



Extending the baseline: Spitzer Mid-Infrared Photometry 
of Globular Cluster Systems in the Centaurus A and 
Sombrero Galaxies 



Lee R. Spitler^* Duncan A. Forbes^ and Michael A. Beasley^ 

^ Centre for Astrophysics & Supercomputing, Swinburne University, Hawthorn, VIC 3122, Australia 
^ Instituto de Astrofisica de Canarias, Via Lactea, E-38200 La Laguna, Tenerife, Spain 



26 June 2008 



ABSTRACT 

Spitzer IRAC mid-infrared photometry is presented for the gfobular cluster (GC) 
systems of the NGC 5128 ("Centaurus A") and NGC 4594 ("Sombrero") galaxies. 
Existing optical photometric and spectroscopic are combined with this new data in 
a comprehensive optical to mid-IR colour catalogue of 260 GCs. Empirical colour- 
metallicity relationships are derived for all optical to mid-IR colour combinations. 

These colours prove to be very effective quantities to test the photometric pre- 
dictions of simple stellar population (SSP) models. In general, four SSP models show 
larger discrepancies between each other and the data at bluer wavelengths, especially 
at high metallicities. Such differences become very important when attempting to use 
colour-colour model predictions to constrain the ages of stellar populations. Further- 
more, the age-substructure determi ned from colou r -colou r diagrams and 91 NGC 5128 
GCs with spectroscopic ages from iBeaslev et all (|2008[) are inconsistent, suggesting 
any apparent GC system age-substructure implied by a colour-colour analysis must 
be verified independently. 

Unlike blue wavebands, certain optical to mid-IR colours are insensitive to the flux 
from hot horizontal branch stars and thus provide an excellent metallicity proxy. The 
NGC 5128 GC system shows strong bimodality in the optical R-band to mid-IR colour 
distributions, hence proving it is bimodal in metallicity. In this new colour space, a 
colour-magnitude trend, a "blue tilt", is found in the NGC 5128 metal-poor GC data. 
The NGC 5128 young GCs do not contribute to this trend. In the NGC 4594 GC 
system, a population of abnormally massive GCs at intermediate metallicities show 
bluer optical to optical colours for their optical to mid-IR colours, suggesting they 
contain extended horizontal branches and/or are younger than typical GCs. Analysis 
of optical to mid-IR colours for a ultra-compact dwarf galaxy suggests its metallicity 
is just below solar. 

Key words: galaxies: spiral - galaxies: star clusters - galaxies: individual (M104, 
NGC 4594, NGC 5128) 



1 INTRODUCTION 

Massive star cluster formation takes place during major star 
bursts. The age substructure of a particular star cluster sys- 
tem therefore provides a historical record of significant star 
formation events in its host galaxy. Globular clusters (GCs) 
are a class of massive star clusters with a mean mass of a 
few IQ^Mq and are notable for their compactness. Observa- 
tional evidences generally suggests the bulk of GC formation 
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took place within roughly 2 gigayears of the Big Bang (see 
references in Brodie & Strader 2006). 

Systems of GCs are known to host a blue and red GC 
subpopulations in terms of their photometric colour (e.g. 
Peng et al. 2006; Strader et al. 2006). This has been in- 
terpreted as evidence for a bimodal metallicity distribution, 
hence most GC systems contain both a metal-poor (or blue) 
and metal-rich (red) GC subpopulation. Given the close 
proximity of the GC formation epoch to the time when the 
universe was free of metals, the existence of two GC metal- 
licity subpopulations likely implies that two distinct phases 
of GC formation took place at different epochs. Obtaining 
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absolute age dates of the two star formation events (and sub- 
sequent ones) can help constrain the detailed formation and 
assembly history of a galaxy. This has encouraged efforts to 
age-date GCs and consequently a number of techniques have 
been employed to ov ercome the age-metallicity degeneracy 
(e.g. IWorthevl[l993 ). 

Analysis of absorption line indices from spectroscopy is 
one such technique and most of the few hundred or so GCs 
with good data are found t o be as old as Galactic G Cs (see 
summary of observations in iBrodie fc StradeJ2006l ). Unfor- 
tunately, the lengthy telescope time required for such analy- 
sis is costly, hence the number of GCs studied in this manner 
remains remarkably low. A n otable exception is the recent 
work of iBeaslev et al.l l|2008l ). who have obtained spectro- 
scopic ages and metallicities for ~ 150 GCs in the massive el- 
liptical galaxy, NGC 5128. Existing spectroscopic studies are 
also biased towards the brightest GCs, an important caveat 
that should be re-empha sised given the discovery of a GC 
colour-magnitude trend i Harris et al. Il2006bl :[ Strader et al.l 
l2006l : ISpitler et al.ll2006l : 



Mieske et al. 



20061 ). 



On the other 

hand, the well-studied Milky Way GC system exhibits no 
such trend with GC luminosity. 

GC colour data can also help overcome the age- 
metallicity degeneracy. This involves a comparison between 
observational colour-colour data for individual GCs and 
theoretical stellar population model predictions. Optical 
to infrared (IR) colours are essential for this work, as 
the degeneracy among optical colours is severe for ages 
> 2Gyr (e.g. Lee et al. 2007a). From analysis of opti- 
cal and near-IR photometry, GC system age substructur e 
has been re ported in NGC 1316 (iGoudfrooii et al.l [200ll ). 
N GC 4365 JPuzia et all I2OO 2I: iKun du et all l2005l: how ever 
c.f. lLarsen, Brodie. fc Strade r 2005; Brodie et al.ll2005h and 
NGC 5846 (jHempel et al.ll2003i ). In most instances, a sub- 
population of "intermediate-age" GCs (~ 2 — 8 Gyr) that 
make up a small fraction (~ 15%) of the total GC sys- 
tem is claimed. A notable exception is the GC system of 
NGC 4365, which might contain ~ 40 — 80% of these younger 
GCs (Puzia et al. 2002; though cf. Larsen et al. 2005 & 
Brodie et al. 2005). The systems for which no age sub- 
structure was reported using near-IR photometry include: 
M87, NGC 4478 (Kissler-Pat ig, Brodie & Minniti 2002), 
NGC 3115 (|Puzia et al.ll2002l ). NGC 4594, NGC 3585 and 
NGC 4472 (Hempel et al. 2007). However, the data samples 
for each of these systems are typically small: on order of a 
few tens of GCs per system. 

Optical to IR colour data is also well-suited for char- 
acterising GC metallicity distributions, hence can help con- 
strain the chemical properties of galaxies at early epochs. 
Such colours span very large dynamical ranges and allow 
more precise metallicity measurements. For example, 0.1 dex 
in metallicity corresponds to ~ 0.03 mag in B-R while 0.1 
dex is ~ 0.07 mag in B-K (according to the preliminary 
models of Chariot & Bruzual 2008 priv. comm.). The IR is 
also insensitive to hot horizontal branch stars. This is a im- 
portant point considering the cautionary work of Yoon, Yi 
& Lee (2006) who suggested GC colour bimodality might 
entirely result from a rapid transition between blue and red 
horizontal branches at intermediate GC metallicities. 

The above issues are addressed here wit h an observa- 
tiona l study of the Spitzer Spac e Telescope (I Werner et al.l 
I2OO4I ) IR Array Camera (IRAC; iFazio et al.ll2004l ) mid-IR 



properties of GC systems in the NGC 4594 and NGC 5128 
galaxies. NGC 4594, the "Sombrero galaxy" is classified as 
an Sa-type spiral, but has a gala xy stellar mass-no rmalised 
number of blue GCs {Tuue; ,Zepf fc Ashmanll 19931 ) compa- 
rable to massive, early-type gal axies (Spitler et al. 2 008). 
With an estimated 1900 GCs (|Rhode fc Zepj |2004| ). the 
NGC 4594 galaxy perhaps hosts the largest GC system 
within the local volume. Recent Hubble Space Telescope 
(HST) Advanced Camera for Surveys (ACS) imaging re- 
vealed bi-modality in the optical colours of the GC system 
(Spitler et al."2006'),^ confirming p r evious work in the optical: 
Gcbhardt & Kisslcr -Pati3 (| 19991 1 . iLarsen. Forbes, fc Brodid 
(|2001i ). and Rhode & Zepf (2004) . Co-added spectra from 14 
N GC 4594 GCs were found to be consistent with old ages 
bv lLarsen et al.l l|2002l) . Comparing V-I and V-K photomet- 
ric data of 26 NGC 4594 GCs to stellar population models, 
Hempel et al. (2007) concluded these GCs are predominantly 
old. 

NGC 5128, "Centaurus A", is a peculiar S O galaxy 
whose GC system shows optical colou r bimodality l|Reikubal 
I2OOII : IPeng. Ford, fc Freemanll2004bl ). Harris et al. (2006b) 
estimated the GC system contains 1550 ± 390 GCs. A 
intermediate- ffl^erf subpopulation making up no more than 
15% of the GC sy stem total was recently reported by 
iBeaslev et al.l l|2008l ). In addition to providing spectroscopic 
eviden ce for the me t al-po or and metal-rich GC subpopula- 
tions. iBeaslev et al.l (|2008l ) also found a small, intermediate- 
metallicity GC subpopulation. 

In the following, the distances to NGC 4594 and 
NGC 5128 galaxies are taken to be m-M= 29.77 ± 0.03 
or 9.0 Mpc Jspitler et al."2006') and m-M= 27.74 ± 0.14 or 
3.5 Mpc HFerrarese et al.ii2007i ). respectively. 



2 DATA 

2.1 Spitzer Data Reduction 

Figure [T] presents IRAC images of the NGC 4594 galaxy 
produced by the Spitzer Infrared Nearby Gala xies Survey 
(SINGS) legacy program l|Kennicutt et al.ll2003D . The IRAC 
third ([5.6]) and fourth ([8.0]) channel images (with filters 
centred on 5.6 and 8.0 microns, respectively) have an over- 
whelming background level and therefore GCs are not ap- 
parent in these images. In contrast, IRAC channels one ([3.6] 
) and two ([4.5]) show many sources over the entire field of 
view. NGC 4594 observations were designed so that the cen- 
tral mosaic at least had 8 x 30s overlapping individual expo- 
sures and the [3.6] field of view is ~ 13' x 25' (~ 34 kpcx65 
kpc). The [3.6] images were dithered by the SINGS group to 
increase the image resolution and astrometric precision. The 
resultant pixel scale of 0.75"pixel~^ is still relatively large 
for this GC system, thus object crowding is problematic. 

The NGC 5128 IRAC image set is associated with the 
Spitzer program #101 (P.I. C. Lawrence) and was prepared 
for analysis with the standard Spitzer archive data pipeline. 
The [3.6] mosaic covers a 27' x 24' (- 29 kpcx24 kpc) 
region slightly ofi'set from the galaxy centre and is shown 
in Figure [2] The mosaic average exposure time is approx- 
imately 72s (some locations on the mosaic have lower or 
higher coverage from the dithering pattern) and the pixel 
scale is 1.2"pixel~^. 
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Figure 1. Spitzor IRAC [3.6] -band image of the NGC 4594 galaxy. GCs candidates from HST optical data (Spitler et al. 2006) that 
had counterparts on the IRAC images are highlighted with yellow circles. 




Figure 2. Spitzer IRAC colour images of the NGC 5128 galaxy. Circle symbols represent the location of confirmed GCs from the 
catalogue of Woodley et al. (2007) that were successfully matched to objects in the Spitzer [3.6] imaging. Red, c yan and yellow circle s 
represent GCs with spectroscopic old, intermediate to young and no ages, respectively according to the analysis of lBeaslev et aL 1 Hooi). 
The lower zoomed inset shows residuals from PSF fits to the [3.6] image presented in the upper inset. 



2.2 Spitzer Photometry 

The software tools used to extract IRAC photometry are 
packaged in the multi -purpose image analysis sof tware ap- 
plication (MOPEX; iMakovoz fc Marleaul l2005h and the 
specific method employed was that recommended in the 
MOPEX documentation. Object detection and photometry 
was carried out on the [3.6] and [4.5] images after the diffuse 
galaxy light was removed by the subtraction of a median- 
filter image. A local object detection threshold of 4(t was 
used. To address object crowding, the apparent magnitude 
of each object was determined by fitting a point spread func- 



tion (PSF) to its profile. An empirical PSF for each IRAC 
dataset was constructed with MOPEX using bright, iso- 
lated objects. The PSF full-width at half-maximums are 1.7 
(2.0') and 2.1 (1.6') pixels for the NGC 5128 and NGC 4594 
[3.6] images, respectively. As demonstrated in Figure [2] the 
PSF fitting residuals are minimal, suggesting the photomet- 
ric method is robust. Fluxes from individual objects were 
calibrated to the Vega magnitude system using the photo- 
metric zeropoints from lReach et al] (|2005l ). Apparent mag- 
nitudes were not corrected for Galactic extinction because 
dust has a negligible impact in this wavelength regime {Al ~ 
A[3,6, < 0.02 mag). 
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Figure 3. Profile fitting completeness results for NGC 5128 (left) 
and NGC 4594 (right). For each panel, the left hand y-axis is the 
[3.6] -band uncertainty and the right hand y-axis is the complete- 
ness fraction. The residuals from the tests (i.e. the absolute value 
of the difference between the input and output magnitudes) are 
shown as the shaded region with the la residual level marked as 
a solid line. Green, solid circles represent objects analysed further 
(A([3.6] )< 0.2 mag), while the open circles were not used. The 
dashed line corresponds to the completeness fraction estimated 
from the tests. 



While the majority of the objects detected in the [3.6] 
-band images are also found in the [4.5] images, the magni- 
tude difference between the [3.6] and [4.5] bands are essen - 
tially negUgible for old GCs (e.g. iBruzual fc Charlotll2003l l. 
For such stellar populations, if the [3.6] and [4.5] photomet- 
ric data had shown comparable measurement uncertainties, 
the two bands might be combined to increase the accuracy of 
their mid-IR photometric measure. In the present case, the 
formal [4.5]-band photometric uncertainties are twice those 
of the [3.6] -band, thus the subsequent analysis only makes 
use of the [3.6] -band data. 

To test the robustness of the Spitzer photometry extrac- 
tion method, images containing artificial objects of known 
magnitudes were produced and their magnitudes were ex- 
tracted using the same method described above. In detail, 
for each 0.2 mag. interval over the observed magnitude range 
of the GCs, an artificial field from the original images field 
was populated with 4000 artifical objects. The same detec- 
tion and PSF-fitting method were applied to these artifi- 
cal fields and the results are summarised in Figure O From 
the residuals between the input and measured magnitudes, 
it is apparent the formal photometric errors estimated by 
MOPEX were slightly underestimated. The formal photo- 
metric errors are therefore taken as the Icr level of the test 
residuals. 



2.3 Optical Photometry 

NGC 4594 GC optical photometry was extr acted from a 
HST A CS image mosaic (B,V, & R bands) bv lSpitler et all 
^2Q0&). The ACS imaging covers a ~ 10' x 7' (~ 26 kpcx 18 
kpc) region, centred on NGC 4594. Optical photometry for 
the NGC 5128 GCs is from the Woodley et al. (2007) cata- 
logue, which is largely drawn from the work of Peng et al. 
(2004). A subset of the NGC 5128 GCs have spectrosc opic 
ages and metallicities measured bv lBeaslev et al.l l|2008l ). In 
both optical datasets, correct ions for Galactic extinction are 
from the DIRBE dust maps (|Schlegel. Finkbeiner. fc Davis! 
Il998l) . 

NGC 4594 GC photometry was derived in lSpitler et al.1 

('2006') using a standard aperture photometry technique. 
They applied a single aperture correction to account for 
light falling beyond the extraction aperture (5 ACS pixels 
or 0.25") and within a radius of 10 ACS pixels (0.50"). A 
constant aperture correction for the entire dataset may not 
accurately represent the GC photometry. More specifically, 
a trend is observed where the more extended objects show 
~ 0.2 mag. more flux than the most compact GC in the 
5 to 10 pixel annulus. The effect is roughly the same in 
each of the B, V and R bands thus the effect was negligi- 
ble in the Spitler et al. (2006) work because it will cancel 
when optical colours are produced. The cancelling occurs 
because the pixel scales of the ACS images are identical and 
filter-dependent PSF variations are small. However, when 
ACS and Spitzer photometry are combined, the aperture 
corrections applied must be performed properly in an abso- 
lute sense to avoid introducing scatter in the ACS-Spitzer 
colours. Such scatter results from the different pixel scales 
of the instruments: GCs are resolved in the ACS images, 
but not in the IRAC images. To mitigate this effect, the 
NGC 4594 optical photometry was re-derived with an ex- 
traction aperture of 10 pixels and adopted throughout the 
following. 



2.4 Optical/Spitzer Cross-matching 

Cross-matching between the optical GC catalogues and 
Spitzer objects was carried out using a matching threshold 
of 1.08" in both the NGC 4594 and the NGC 5128 images. In 
the NGC 4594 dataset, to decrease the chances of a nearby 
(projected) object effecting Spitzer flux measurements, all 
GC candidates with at least one nearby (< 1.6") object in 
the ACS images of a comparable magnitude (within +2.5 
B-band magnitudes of the GC candidate) were culled from 
further analysis. A total of 19 NGC 4594 GC candidates 
were not considered for this reason. NGC 5128 GCs were 
not subject to this procedure because optical imaging is not 
publicly available. A single NGC 4594 GC candidate shows a 
signiflcant excesses of [3.6] light expected for GCs (Spitler et 
al. [2006] ID#538; with V = 22.98 and [3.6] = 17.81). This 
GC candidate was removed from further analysis as a pos- 
sible contaminate, although flux confusion from an object 
undetected on the ACS imaging (e.g. a distant red galaxy) 
cannot be ruled out. 

A photometric error selection of A([3.6] ) ^ 0.2 mag 



was applied. This roughly corresponds to Mrj 



[3.6] 



^10.2 



and M 



[3.6] 



-11.4 mag. A smaU fraction of the NGC 5128 
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Figure 4. Colour-magnitude diagrams for NGC 5128 (top) and 
NGC 4594 (bottom) GCs. A UCD associated witli tlie NGC 4594 
galaxy is given in the two lower diagrams at bright magnitudes. 
Open circles in the two left plots have no corresponding match on 
the Spitzer imaging. For NGC 5128, red, blue and black circle s 
represent GCs with spectroscopic ages from lBeaslev et al.1 (I2OO8I ') 
with ^ 8Gyr, < 8Gyr and unknown ages, respectively. GC colours 
show bimodal distributions. 



Figure 5. Same layout as Figurc|4l with the corresponding colour 
histograms given. Dashed red histograms show NGC 5128 GCs 
with old 8Gyr) ages and the green histogram in the lower left 
plot shows the B— R colour histogram of NGC 4594 GCs matched 
with Spitzer photometry. Colour bimodality is most apparent in 
the optical, NGC 5128 B-[3.6] and NGC 5128 R-[3.6] distribu- 
tions. 



optical photometry is very uncertain, therefore an additional 
selection of A(B— R ) ^ 0.1 mag. was applied. 

After this selection, 76 and 184 GC candidates with 
both optical and Spitzer photometry remained in the 
NGC 4594 and NGC 5128 datasets, respectively. Of the 
matched NGC 5128 GCs with good spec troscopic age mea- 
surements (91) from lBeaslev et al.l l|2008l ). 11 or ~ 12% show 
ages less than 8 Gyr. All [3.6] photometry is available from 
the first author and is posted on the CDS. 

A very luminous, compact object (V — 17.46) that is 
well-resolved on the HST ACS imaging, was also detected 
on the [3.6]-band imaging and has [3.6] = 14.19. A radial 
velocity measurement from the spectroscopic sample of Hau 
et al. (2008, in prep.) suggests this object is associated with 
the NGC 4594 galaxy, hence is considered a ultra-compact 
dwarf galaxy (UCD) candidate. It is included in the follow- 
ing analysis when appropriate. 



3 RESULTS/DISCUSSION 

3.1 Optical/Mid-IR Colour Distributions 

3.1.1 Bimodality 

GC colour-magnitude diagrams are presented in Figure |4] for 
both GC systems. Colour histograms are given in Figure [S] 
Bimodality in the optical to Spitzer colours is clearly evident 
in the Figures and the colour "gap" between the GC sub- 
populations is located at B— [3.6] 3.4 and R— [3.6] ~ 2.1 
mags. 

The KMM algorithm of lAshmanT Bird, fc Zepj (|1994| ) 



was used to quantify how likely these datasets are drawn 
from a bimodal distribution instead of a unimodal one. 
When run on the NGC 5128 B-[3.6] data, the KMM output 
strongly supports colour bimodality at a highly significantly 
level (see Table [l| . If a KMM analysis is performed only 
on GCs that are spectroscopically determined to be old, the 
statistical confidence for B — [3.6] bimodality remains strong. 

The KMM results support the idea that the NGC 4594 
B— [3.6] distribution is well-represented as two Gaussian 
compared to a single Gaussian. However, the KMM B— [3.6] 
colour peaks are ~ 0.5 mag redder than what can be in- 
ferred from the colour-magnitude diagram in Figure |4] Two 
features found by Spitler et al. (2006) in the NGC 4594 
colour-magnitude parameter space likely influence the KMM 
results: the existence of bright GCs with colours consistent 
with the subpopulation gap and blue GCs show a colour- 
magnitude trend (the "blue tilt"; see H3.1.2p . thus brighter 
blue GCs will tend to "fill-in" any gap between two colour 
subpopulations. A more complete sampling of the GC lumi- 
nosity function is required before the B — [3.6] KMM output 
is to be trusted. 

This analysis confirms GC colour bimodality previously 
found using optical photometry (see references in the Intro- 
duction). It is apparent that the larger baseline of the opti- 
cal to Spitzer colours helps identify the GC subpopulations 
more clearly than when only optical photometry is used. 

3.1.2 The Blue Tilt 

Spitler et al. (2006) discovered a trend among the blue 
NGC 4594 GCs, where the colours become redder with in- 
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Table 1. Results from a KMM analysis of the GC B-R , B-[3.6] 
, and R— [3.6] distributions. The p-stat value is the bimodality 
gauge produced by the KMM algorithm (see text). The last two 
columns are the mean colour of each GC subpopulation estimated 
by KMM. To prevent data outliers from influencing the KMM 
tests, the NGC 4594 colours were restricted to B— [3.6] ^ 4.8 and 
the NGC 5128 data to 2.4 < B-[3.6] ^ 4.2. The R-[3.6] data 
is bounded by 1.3 < R— [3.6] < 3.2. The heteroscedastic version 
of the KMM test is employed for the B— [3.6] and R— [3.6] colour 
distributions. 
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Figure 6. Spitzer to optical colour-magnitude diagrams. Blue 
and red points correspond to objects with B-R ^ 1.3 and B-R 
> 1.3, respectively. This split on B-R corresponds to the GC 
subpopulation division used by Spitler et al. (2006). Lines fitted 
to the subpopulation colour-magnitude data are given. A blue GC 
colour-magnitude trend (a "blue tilt") is found in the NGC 5128 
data, while incomplete luminosity coverage of the NGC 4594 data 
limits such analysis. 



Table 2. Linear fits to GC subpopulations. The NGC 4594 data 
wore restricted to B-[3.6] ^ 4.8 and the NGC 5128 to 2.4 < 
B— [3.6] ^ 4.3 before the fits were performed. 
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red 
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68 
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old 


red 
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old 
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red 


46 


-0.293 ±0.020 


9.06 ±0.34 


A[3.6] ^ 0.2 


blue 


28 


-0.277 ±0.025 


7.89 ± 0.42 



creasing GC luminosity. Harris et al. (2006a) pointed out a 
similar trend is apparent in a figure of Harris et al. (2004) 
containing NGC 5128 photometry. To explore this so-called 
"blue tilt" (Strader et al. 2006) in the optical to Spitzer pa- 
rameter space, the GC systems was divided into a blue and 
red subpopulation at B— [3.6] ^ 3.4 and B— [3.6] > 3.4, re- 
spectively. Each subpopulation was then fitted with a line to 
the B— [3.6] versus [3.6] information (weighted by the formal 
photometric uncertainties). Parameters are in Tableland 
the fits themselves shown in Figure [6l 

The linear fits to the NGC 5128 GC subpopulation data 
confirm the qualitative assessment of existing optical pho- 
tometry by Harris et al. (2006a) and suggest the blue sub- 
population shows a very significant (~ 9cr) B— [3.6] versus 
[3.6] colour-magnitude gradient, while the red subpopula- 
tion shows no such trend. Fitting only to old NGC 5128 
GCs yields a similar outcome. Converting the NGC 5128 
to a metallicity-mass proportionality with the empirical 
colour- metallicity transformations described in Sec. 13.21 the 
NGC 5128 blue GCs show Z oc M°-^^ The strength of this 
effect is among the lowest reported thus far. The approxi- 
mate mass range covered by the Spitzer data is ~ 5 x 10^ 
to IO'^Mq. 

Fits to the NGC 4594 [3.6] versus B-[3.6] data suggests 
both GC subpopulations show a formally significant (~ IOct) 
colour- magnitude trend (see the fits in Fig. IS}. A red sub- 
population trend was not found in the optical analysis by 
Spitler et al. (2006). A better sampling of the GC luminos- 
ity function (which currently only covers ~ 10*^ to lO'' Mq) 
is likely required before such trends are constrained with 
Spitzer mid-IR photometry. 

The cause of blue tilt is not currently understood 
though a number of ideas have been proposed (see discus- 
sion in Strader et al. 2006; Harris et al. 2006a; Mieske et al. 
2006; Spitler et al. 2006). The present work rules out age 
effects as the cause of the NGC 5128 blue tilt. If the trend is 
due to increasing metallicity with luminosity, an analogous 
trend for red GCs should be a factor ~ 10 times shallower 
than the gradient observed among the blue GC subpopu- 
lation (Spitler et al. 2006), hence a good reason for it not 
being observed if it is indeed present. Determining whether 
or not a metallicity-induced red tilt actually exists would 
clearly be aided with the long baseline offered by Spitzer to 
optical colours. 



Spitzer IRAC Photometry of GC Systems 7 




0.4 0.6 0.8 1 l.aO.5 1 1.5 1 1.5 2 2.33.2 0.4 0.6 




Figure 7. Colour- metallicity diagrams of NGC 5128 G Cs with ages > 8 Gyr and metallicity uncertainties of < 0.15 dex. NGC 5128 
metallicities ([M/H]; scale is from Harris 1996) are from iBeaslev et al.l 1 12008 ). Black lines are the best fitted linear relationships to the 
data using a least-squares fit (clipping outliers, as appropriate). The parameters of these fits are presented in Table l3l Predictions from 
five stellar population models are shown. 



3.2 Optical-Spitzer/Metallicity Relationships 

Figure [7] provides all Johnson/Cousins optical to optical and 
optical to Spitzer colour permutations against GC metal- 
licity for the ^ 8 Gyr NGC 5128 GC data. The colour- 
metallicity data are well-represented with linear fits, which 
are presented in the same Figure and detailed in Table [S] 
These relationships are only valid for stellar population 
with ages > 8 Gyr and metallicities of —1.8 < [M/H] 
< -1-0.0. Individual NGC 5128 GC metalhcities were de- 



rived in lBeaslev et al.l (|2008l ) and the t echnique employed is 
briefly summarised here. For each GC.'Be aslev et al] (|2008l ) 
took every Lick index measured and transformed it into a 
metallicity value using empirical relationships derived from 
Galactic GC data. The Tukey biweight of the all Lick index 
metallicity values is taken as the final metallicity ([M/H]) for 
a particular GC. The [M/H] scale is that of Harris (1996) 
for the Galactic GC system. Individual NGC 5128 ages were 
estimated with fits to the Lee & Worthey (2005) SSP mod- 
els. 
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Table 3. Empirical NGC 5128 GC optical to mid-IR colour- 
mctallicity relations. Valid for use with old (> 8 Gyr) single stel- 
lar populations spanning a metallicity range of —1.8 < [^'I/H] 
< +0.0. Metallicity scale is that of the Harris (1996) Galactic GC 
catalogue. Optical colours are corrected for Galactic dust extinc- 
tion using the maps of Schlegel et al. (1998). 
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The optical to [3.6] colours show a stronger correlation 
with metallicity compared to optical-optical colours (see for- 
mal significance of the fitted slopes in Table |3} . To what 
degree this reflects the possibility that a linear fit is more 
appropriate for optical to [3.6] colours or else the fact that 
[3.6] is relatively insensitive to second parameters (such as 
age and horizontal branch morphology) is unknown. Cer- 
tainly the latter point has some influence and the larger 
dynamical range provided by optical to [3.6] colours makes 
them a preferred metallicity proxy over optical to optical 
colours. 

Certain SSP models (e.g. Bruzual & Chariot 2003 here- 
after BC03; Chariot & Bruzual 2008 priv. comm., CB08; 
Maraston 2005, M05; S. Yoon 2006 priv. comm. , YQ6) pro- 
vide optical to Spitzer colour predictions for a SSP of a given 
age and metallicity. As GCs are single stellar populations 
whose stars show a single age and metallicity, they pro- 
vide an ideal metric to contrast different model predictions 
for old stellar populations (e.g. Proctor, Forbes, & Beasley 
2004; M05; Mendel et al. 2007; CantieUo & Blakeslee 2007; 
Cohen et al. 2007; Pessev et al. 2008). Figure [7] provides 
the 12 Gyr colour-metallicity prediction for five SSP mod- 
els: BC03, CB08 (preliminary update to BC03), M05 (blue 
horizontal branch), Raimondo et al. (2005; SPoT), and Y06. 
Note for the SPoT models, the L-band predictions are used 
in place of the [3.6]-band, which are currently not available. 
For Y06, 13 Gyr models are used to best illustrate their 
predicted non-linear structure due to changing horizontal 
branch (HB) morphology (see Sec. I3.2.T|) . 

It is apparent in Figure [7] that all model colour pre- 
dictions are roughly consistent with each other. When the 



models are compared to the observed NGC 5128 GC colours, 
it is apparent at [M/H] < —1.5 the models predict steeper 
[M/H]-colour gradients. This is especially pronounced in 
colours containing redder optical wavelengths such as R—, 
/— , and [3.6]-band. The apparent internal consistency be- 
tween the models, whose treatment of HB stars differs (e.g. 
blue vs. red HBs), argues against HB morphology as the 
primary cause for the discrepancy between the data and 
models. The BC03 and CB08 models tend to predict redder 
optical-[3.6] colours for a given metallicity at [M/H] > —0.5. 

The empirical linear colour-metallicity relations pro- 
vided in Table 13] should prove useful when the observed pho- 
tometric uncertainties are comparable to the uncertainties 
of the NGC 5128 GC photometry. An attempt to better 
constrain the SSP models with the observations is made in 
Section [3]3TT] where the R— [3.6] colour is used as a model- 
independent proxy for metallicity. 

3.2.1 Does GC Colour Bimodality Imply Metallicity 
Bimodality? 

Yoon et al. (2006; see also Richtler 2006) recently demon- 
strated that the interpretation of a bimodal optical colour 
distribution as evidence for two subpopulations of GCs de- 
pends critically on the relationship employed to translate 
colour into metallicity. According to their SSP models (Y06), 
the transition from a blue to red horizontal branch (HB) at 
intermediate metallicities causes a change in slope of this 
relation. This change in slope means that an observed bi- 
modal colour distribution may in reality be resultant from 
a single, unimodal metallicity distribution. 

To determine how much of an effect HB stars have on 
the observed GC colour bimodality, colours derived from 
wavebands that are insensitive to HB stars (e.g. IR) can be 
compared to those colours that are (e.g. B-band). Indeed, 
simulations by CantieUo & Blakeslee (2007) of the SPoT 
models show that certain optical/near-IR colours (e.g. V-H 
and V-K) are preferred for circumventing any biases caused 
by any non-linear structure in the colour-metallicity rela- 
tionship. They also demonstrate that the non-linear struc- 
ture is not consistent between the models they examined 
(they compared 5 different SSP models), thus any effect like 
that described in Yoon et al. (2006) is strongly dependant 
on the particular colour-metallicity relationship adopted. 

S. Yoon (2006, priv. com.) kindly provided [M/H] pre- 
dictions from the same SSP models used in Yoon et al. 
(2006) for three colours (B-R , B-[3.6] , and R-[3.6] ), 
which are shown in Figure [T] The inflection point from a 
changing horizontal branch type at [M/H] ~ —0.8, is appar- 
ent in their B— R model predictions. The predicted [M/H] 
versus R— [3.6] relationship of Y06 in Figure [7] is approx- 
imately linear, suggesting horizontal branch stars are no 
longer playing an important role at these wavelengths. S. 
Yoon (2006, priv. com.) included a cautionary note that, in 
general, mid-IR stellar libraries are poorly populated, thus 
any model prediction at these wavelength are not well con- 
strained. Nevertheless, the relative insensitivity of R— [3.6] 
to HB stars predicted by Y06 makes them a useful colour 
combination to test the scenario of Yoon et al. (2006). 

Analysis in Sec. 15X1] of the old NGC 5128 GC R-[3.6] 
data provides strong evidence for colour bimodality. Since 
this colour is insensitive to horizontal branch stars, these 
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observations suggest colour bimodality is likely not caused 
by a variable horizontal branch morphologies, as suggested 
by Yoon et al. (2006). Given the old ages of these GCs, 
the NGC 5128 colour distribution implies a bimodal metal- 
licity distribu t ion. T his is in agreement with the results of 
iBeaslev et al.1 (|2008h using spectroscopic data of NGC 5128 
GCs. 

Although the NGC 4594 optical to mid-IR colour data 
are not strongly bimodal, there is a good reason why this 
is expected. In Figure |4l it is apparent that the NGC 4594 
GCs with Spitzer photometry only cover ~ 2 magnitudes of 
the GC luminosity function. As the brightest, blue GCs are 
biased towards redder colours, they effectively "fill-in" the 
subpopulation gap, hence making any inherent bimodality 
more difficult to detect. The NGC 5128 data samples almost 
twice the magnitude range, hence the blurring of the sub- 
population data at high luminosities is less critical for this 
analysis. 

Kundu & Zepf (2007) concluded the observed bimodal 
I-H colour distribution of very bright GCs in the giant el- 
liptical M87 as evidence for a bimodal metallicity distri- 
bution with a similar argument to that used here. Strader, 
Beasley & Brodie (2007; also Beasley et al. 2008) tackled the 
GC metallicity bimodality question without the use of the- 
oretical models and employed empirical relationships they 
derived between metallicity and absorption line indices of 
Galactic GCs. Although their technique did not explicitly 
account for variable HB morphology effects, they provide 
compelling evidence that luminous GCs in the giant ellipti- 
cal NGC 4472 are drawn from a bimodal metallicity distri- 
bution. 

The optical to Spitzer colours produced here provide 
strong evidence for GC system metallicity bimodality in 
NGC 5128, and support recent results for bimodal GC 
metallicity distributions in certain giant elliptical galaxies. 



3.3 Colour-Colour Diagrams 

3.3.1 SSP Model Constraints 

The use of optical to Spitzer colours to test the colour predic- 
tions from different SSP models has certain advantages over 
metallicity. The term "metallicity" is an ill-defined quan- 
tity and can vary between different SSP models (see discus- 
sion in Beasley et al. 2008) . Colours, on the other hand, are 
straight-forward to compute and therefore should be directly 
comparable between SSP models. Optical to Spitzer colours 
are very good metallicity proxies and provide a higher pre- 
cision due to increased dynamical ranges ( i]3.2p . hence they 
provide an excellent quantity for comparing theoretical pre- 
dictions. 

Figure |8] shows the same colours of Figure [7] instead 
plotted against the R— [3.6] colour. The R— [3.6] was chosen 
over e.g. U-[3.6], B-[3.6], etc., for the following reasons: (1) 
the NGC 4594 data can be included in the analysis; (2) 
it is relevant to the age-substructure analysis of i|3.3.2l (3) 
R-[3.6] shows most significant statistical correlations with 
[M/H] (see Table |3]). The following analysis holds true if 
B— [3.6] is used instead as the basis. Model predictions for 12 
Gyr ages are presented in Figure [S] The mean spectroscopic 
age of the NGC 5128 blue and red GC subpopulations are 
11.1 and 10.8 Gyr, respectively, each with a la age spread 



of ~ 0.9 Gyr (|Beaslev et al.ll2008l ). Analysis of roughly two 
dozen NGC 4594 GCs suggests they have old ages (Larsen 
et al. 2002; Hempel et al. 2007). The following analysis does 
not change if the NGC 4594 data is not used. 

Overall the agreement between the model predictions 
and the data is very good. In Figure O it is apparent that 
GC data generally reach lower values of R— [3.6] than pre- 
dicted by the models, perhaps indicating the predicted [3.6] 
flux is underestimated at low metallicities. To better facili- 
tate a comparison, the data and model residuals from the fits 
to the data (shown in Figure|8| are plotted in Figure[9l Visu- 
ally inspecting each optical wave-band in turn, the following 
observational constraints can be deduced from Figure [§1 

• [/—band. The BC03 and CB08 colour predictions best 
reflect the shape of the data. The SPoT and especially M05 
models predict bluer colours (more band flux) at low- 
metallicities and redder colours (less [/—band flux) at high- 
metallicities for colours including this band. CB08 provides 
the best match to the data at low- metallicities. 

• B— band. The i?— band model residuals resemble those 
of 17— band. M05 predicts redder colours (less B— band flux) 
at high metallicities compared to the other models and the 
data. The other models actually show a similar residual 
structure to M05, although it is less pronounced. Again, 
CB08 matches the data at the low metallicities. At inter- 
mediate metallicities the SPoT models match the data well 
and show some resemblance to the Y06 models in the B-R 
diagram. 

• y-band. For this filter, the M05 and SPoT models fair 
best, while BC03 and CB08 models evidently predict less V- 
band fiux compared to the data at all metallicities. This best 
illustrated in the V-R, V-I and V-[3.6] diagrams, but is also 
apparent when comparing the BC03 and CB08 predictions 
in the U-V diagram to U-B. 

• i?— band and 7— band. The predictions and data appear 
to be in good agreement with each other, as demonstrated 
by the R-I and I-[3.6] diagrams. If this is confirmed, it also 
suggests the data and model predictions at [3.6] are fairly 
consistent. 

To conclude, optical to Spitzer colours are shown to 
provide a useful basis to help constrain SSP model pre- 
dictions at optical wavelengths. A comparison between the 
NGC 5128 and NGC 4594 GC data suggests SSP model pre- 
dictions at bluer wavebands {U, B, and V) show the largest 
discrepancies with the GC data compared to the red bands 
{R and /). 



3.3.2 Photometric Ages? 

GC colour-colour plots which include K-band photometry 
(V-I versus V-K is popular for GC work) have been used to 
look for age-substructure in GC systems (see references in 
the Introduction). In theory, this technique works because 
an optical to optical colour (e.g. V-I) is slightly more sen- 
sitive to age, while an optical to K-band colour (e.g. V-K) 
is sensitive to metallicity. Th e NGC 5128 GC spectroscopic 
analysis from iBeaslev et al.l (|2008l ) enables the first large 
comparison between spectroscopic GC ages and those de- 
rived from optical and IR colour-colour diagrams. Although 
the spectroscopic ages are not free from uncertainties, they 
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Figure 8. Colour-colour diagrams of NGC 5128 GCs with ages 8 Gyr (black circles) and NGC 4594 GCs (blue circles). The X-axis 
colour label is given within each plot. Only data with X-axis colour errors of < 0.1 mag are shown. SSP models are as in Figure[7] Lines 
are fits to the GC observations with: C = a(R— [3.6] )+/3, where C is the colour on the X-axis. Mean photometric errors of the GC data 
are illustrated with the crosses in each diagram. 



better break the age-metallicity degeneracy (e.g. Worthey 
1994). 

Li & Han (2008) investigated the relative sensitivi- 
ties of all possible colour-colour pairs among the standard 
BVRIJHKugriz band-passes using the SSP of Bruzual & 
Chariot (2003). They concluded that the colour pairs, R- 
K versus m-R and r— K versus u — r, are most capable of 
partially breaking the age-metallicity degeneracy. The clos- 
est colour combination in the NGC 5128 GC catalogue to 
the preferred pair found in Li & Han (2008) is R— [3.6] ver- 



sus U-R. Analysis in Section 13.3.11 of Figure |9] suggests the 
CB08 provide the best match to old NGC 5128 GCs in this 
colour-colour parameter space. GCs with reliable spectro- 
scopic ages are presented in Figure [TOl with the 12-, 7-, 5- 
and 2-Gyr SSP model predictions from CBQ8. 

It is immediately apparent in Figure [10] that the young 
NGC 5128 GCs determined from spectroscopy show no 
systematic difference to the old GCs. Even the youngest 
GCs (c./. solid blue to open red circles) occupy the same 
parameter-space as the old GCs. Relatively large photo- 
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Figure 9. Data and model residuals from fits to the GC data in Figure |8] Positive residuals correspond to colours redder than the fit 
to the GC data. See text for discussion. 



metric uncertainties and the inherent difficulty in breaking 
the age-metallicity degeneracy using colour-colour diagrams 
provide the most likely explanation f or the discrepa n cy be - 
tween the present findings and the iBeaslev et al.l (|2008l ) 
spectroscopic ages. 

If spectroscopic age information were not available for 
individual GCs, it might be tempting to tag a few of the 
more luminous metal-rich NGC 5128 GCs in Figure [10] as 
candidate young GCs because their colours coincide with a 
younger region of the SSP model grids. Clearly, such a ven- 
ture would prove erroneous once the spectroscopic ages were 
considered, suggesting any apparent age-substructure in 



colour-colour space must be investigated further with spec- 
troscopic data. Furthermore, had another model been em- 
ployed for the colour-colour analysis (without the spectro- 
scopic ages), the interpretation would have varied tremen- 
dously. For example, the region where metal-rich GCs fall in 
Figure [TOl matches with the < 5 Gyr predictions of M05. It 
is therefore apparent that accurately age dating GCs with 
colours alone is limited by the high level of systematic un- 
certainty in SSP models. 

With these caveats in mind, the NGC 4594 GCs are 
presented in R— [3.6] versus B-R colour-colour space in Fig- 
ure 111! Analysis of Figure [9] in Section 13.3.11 provides no 
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Figure 10. NGC 5128 GC R— [3.6] vs. U— R colour-colour diagram. To highlight the most accurate photometry, the circles are inversely 
proportional to photometric uncertainties. Blue and red circles correspond to GCs with spectroscopic ages of < 8 and ^ 8 Gyr, respectively. 
GCs with ages < 4 Gyr are shown as solid circles. CB08 SSP model 2, 5, 7 and 12 Gyr predictions are presented and the arrows indicate 
the general directions in the models for decreasing age and metallicity. A dust extinction arrow of Ay = 0.1 mag. is given on the lower 
right. 



strong preference for a particular SSP to be employed in 
the R— [3.6] versus B-R diagram, thus the CB08 models 
used in Figure [10] are used again. It was also shown in 
Section 13.3.11 that, relative to the data, all models pre- 
dict lower B— band fluxes relative to the GC data at high- 
metallicities, hence GCs will appear younger than they ac- 
tually are on the R— [3.6] versus B-R diagram. Indeed, in 
Figure [TTIthe NGC 5128 GCs with high metallicities tend to 



occupy younger regions in the SSP model grid, even though 
they have spectroscopic ages of ^ 8 Gyr. 

The NGC 4594 GCs show a significant (i.e. greater than 
the photometric uncertainties) spread of colours at R— [3.6] 
> 2.5. This could result from an intrinsic age substructure 
among the metal-rich GCs. Given the possibility that the 
SSP models under-predict the _B— band flux at high metal- 
licities, makes any attempt to interpret the apparent colour- 
spread as an age effect very uncertain. If the data and mod- 
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Figure 11. NGC 4594 and NGC 5128 GCs R— [3.6] versus B— R colour-colour diagram. Layout is like Fig. 1101 and again the circles are 
inversely proportional to photometric uncertainties. NGC 4594 GCs are shown as green circles. Bright ([3.6] < 16.6 mag.) NGC 4594 GCs 
are highlighted in magenta and the UCD associated with NGC 4594 is blue. Although the NGC 4594 GCs occupy the ~ 7 Gyr region 
on the SSP grids, analysis in Section 13.3.11 suggests the model predictions at high metallicities are too red in B— R , thus no significant 
age-substructure can be reported. The Y06 model predictions are presented in yellow for a 13 Gyr (curvy line) and 11 Gyr SSP. See text 
for discussion. 



els were taken at face value, the CB08 predictions suggest 
the metal-rich GCs span an age range of ^ 5 — 12 Gyrs. 
Previous studies have found no evidence for a population 
of intermediate-aged GCs. Larsen et al. (2002) analysed ab- 
sorption line indices of ~ 6 metal-rich NGC 4594 GCs and 
found ages of > 10 Gyr. Hempel et al. (2007) analysed K- 
band data of ~ 1 dozen NGC 4594 metal-rich GCs in a 
similar manner and found no evidence for age substructure. 



At R— [3.6] ~ 2.1 in Figure [TT] a subpopulation of "ab- 
normally luminous" blue GCs (Spitler et al. 2006) are high- 
lighted. In colour-colour space, these GCs show bluer optical 
colours than their fainter GC counterparts, possibly suggest- 
ing they have younger ages of approximately 5 — 7 Gyrs. 
These will be discussed further in the Section [3.3.31 below. 

Shown as a blue circle in Figure[TT]is the UCD candidate 
associated with NGC 4594 (see i|2.4p . Despite the excellent 
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photometric data on this object, the discrepancy between 
the data and SSP models makes it difficult to age-date the 
UCD with just photometry. However, it can be concluded 
that the UCD shows colours consistent with NGC 4594 
metal-rich GCs, suggesting its metallicity is just a few tenths 
of a dex below solar metallicity. 

3.3.3 Hot Horizontal Branch Stars 

From intermediate- to low-metallicities, the NGC 4594 GCs 
in Figure [11] show a non-linear structure in colour-colour 
space that roughly resembles the shape of Y06 SSP 13 
Gyr model predictions (given in the same figure), whose 
inflection-point results from a changing HB morphology (de- 
scribed in Section [323 • Interestingly, the NGC 4594 GCs 
at the point where the YOG inflection point becomes appar- 
ent (B— R ~ 1.3 and R— [3.6] ~ 2.1) are dominated by the 
abnormally luminous blue GCs discussed in Section [3.3.2l If 
these objects were omitted, the resemblance to the Y06 13 
Gyr inflection-point disappears completely. 

Although evidence contrary to the scenario of Yoon et 
al. (2006) has been presented (see Sec. I3.2.1|l . it does not 
rule out their model prediction regarding the effects from 
a changing HB morphology. Testing the YQ6 models with 
observations is difficult because few GCs fall into the gap 
between the GC metallicity subpopulations - exactly where 
the predicted inflection point is predicted by Yoon et al. 
(2006). Furthermore, a strong dependence on age (compare 
the 13 and 11 Gyr YOB predictions in Fie:.lll|) makes it possi- 
ble to tweak models to the observations by carefully picking 
the exact ages of these "old" GCs, since old GCs generally 
have uncertainties of a few Gyrs. Thus the present observa- 
tions can provide no constraint on the possible presence of 
an inflection point as predicted by the Y06 models. 

Another possibility is that these GCs host an exotic HB. 
Certain Galactic GCs host HBs that show a unusually large 
number of very hot stars, giving an extended appearance to 
the HB in individual GC colour-magnitude diagrams. These 
extended HBs tend to be found in more massive GCs in the 
Galaxy (Recio-Blanco et al. 2006; Lee et al. 2007b) and they 
will naturally emit more light at bluer wavelengths com- 
pared to GCs with normal HBs. While the R— [3.6] colours 
of the NGC 4594 luminous blue GCs should not be effected 
by such stars, the B— R colours might be. This combination 
could lead to a situation resembling the data in Figure 1111 

To determine whether B— R colours are sensitive to ex- 
tended HB morphologies. Figure [12] highlights Galactic GCs 
with extended HBs (classified as such by Lee et al. 2007b) . In 
the Figure, it is apparent that GCs with unusually extended 
HBs show colours that are roughly consistent with those 
GCs with normal HBs. This is in agreement with Smith & 
Strader (2007) who found no residual HB trend with optical 
colours (although they didn't analysis B-R colours) after 
correlations with metallicity were removed. 

In Figure [T2]there are two metal-rich Galactic GCs with 
extended HB morphologies (NGC 6388 and NGC 6441; Rich 
et al. 1997; Busso et al. 2007) that have B-R colours (~ 1.3) 
similar to the abnormally luminous NGC 4594 blue GCs. 
The coincidence between these GCs in this parameter space 
and the systematically lower values of B— R for their ob- 
served metallicities (as observed in the NGC 4594 GCs) sug- 
gests the presence of extended HBs in these NGC 4594 GCs 
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Figure 12. A comparison between Galactic GCs and NGC 4594 
GCs. Galactic data is from Harris (1996) and is restricted to GCs 
with reddening values E(B-V) < 0.5. The NGC 4594 metallicities 
were estimated using their R— [3.6] colours and the [M/H] trans- 
formation of Tab. |3] Galactic GCs with strongly extended HB 
morphologies (from Lee et al. 2007b) are presented with solid, 
blue triangles. Moderately extended HBs are highlighted with 
a cyan triangle. NGC 4594 GCs brighter and fainter than [3.6] 
= 16.6 are presented as solid magenta and open green circles, 
respectively. At the B— R colour of the NGC 4594 massive GCs 
(B-R ~ 1.3), two Galactic GCs (NGC 6388 and NGC 6441; see 
Rich et al. 1997 and Busso et al. 2007) with strongly extended 
HBs show systematically bluer colours than other Galactic GCs 
of a similar [M/H]. 



might be remotely plausible. UV observations are needed for 
confirmation (e.g. Sohn et al. 2006). 



4 CONCLUSIONS 

In this work, a catalogue of 260 GC candidates with good 
optical and Spitzer [3.6] photometry is presented, enabling, 
for the first time, an examination of the optical to mid- 
IR properties of a large sample of GCs. Following previ- 
ous examples in the literature using K-band photometry, 
the age-substructure of these GC systems is examined with 
colour-colour diagrams including the [3.6]-band. A spread of 
ages in the NGC 5128 metal-rich GC subpopulation might 
have been reported from photometry alone. However, the 
youngest GC s foun d in the large spectroscopic study of 
iBeaslev et al] l|2008l ) do not occupy a young region in colour- 
colour space and actually show photometric properties that 
are consistent with the old NGC 5128 GCs. This observa- 
tion likely demonstrates that the age-metallicity degeneracy 
is indeed difficult to overcome when using only typical pho- 
tometric observations. Spectroscopic age confirmation is re- 
quired for any age-substructure to be interpreted from old 
stellar population in such colour-colour diagrams. 

For the NGC 4594 GC system, for which no large sample 
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of GCs with spectroscopic ages exists, the metal-rich GCs 
generaUy fall onto the the 5-7 Gyr SSP model predictions 
in B— R versus R— [3.6] . However, an empirical comparison 
between data and SSP models (see below) suggests the mod- 
els likely under-predict the _B— band flux relative to the data 
at these metallicities. Until this discrepancy is understood, 
the apparent subpopulation of intermediate-aged NGC 4594 
metal-rich GCs will remain unconfirmed. 

Another notable feature in the NGC 4594 GC colour- 
colour diagram (Fig. is an apparent non-linearity traver- 
sal of the parameter space, which resemble the SSP model 
predictions of S. Yoon (2006, priv. comm.; Y06). Driving this 
resemblance are a subpopulation of very luminous and ex- 
tended GC candidates first noted by Spitler et al. (2006). 
These objects have B— R colours similar to two Galactic 
GCs that host very extended HB morphologies (Rich et al. 
1997; Busso et al. 2007). These Galactic and the luminous 
NGC 4594 GCs are bluer in B— R than expected for their 
metallicities, perhaps suggesting they share similar HB mor- 
phologies (see Fig. [T2|l. 

Another possibility is that the massive NGC 4594 GCs 
are ~ 5 Gyr younger than the bulk of the GC system. 
This would explain their high luminosities, but not the ex- 
tended sizes they tend to show (Spitler et al. 2006). A dwarf 
nuclei with a young stellar component that accreted onto 
NGC 4594 might match these characteristics, although test- 
ing this scenario with simulations is required. 

Using the excellent mass proxy, [3.6] , and the largest 
colour dynamical range studied thus far for subpopulation 
colour-magnitude trends, the present work demonstrates the 
NGC 5128 blue GC subpopulation shows redder colours at 
higher luminosities. The NGC 4594 sample does not span a 
large enough range of the GC luminosity function to con- 
strain such trends. The NGC 5128 blue GCs show a mass- 
metallicity proportionality oi Z (x M"'^^; noticeably weaker 
compared to other galaxies. This confirms initial suspicions 
of Spitler et al. (2006) that the strength of the blue tilt de- 
creases with the host galaxy luminosity. No such trend is 
found among the red GCs in NGC 5128. Thel Beaslev et al.1 
ifiooa ') NGC 5128 GC spectroscopic ages aUow, for the first 
time, the blue tilt to be studied independent of any young to 
intermediate-aged GC. It is concluded that age substructure 
in not likely playing a role in the blue tilt. 

The scenario of Yoon et al. (2006) is tested with the 
R— [3.6] GC data. Such colours are insensitive to hot HB 
stars (see the Y06 models in Fig [7| , hence can determine 
whether GC colour bimodality implies metallicity bimodal- 
ity without effects from a rapidly changing HB morphology. 
While the NGC 4594 sample is too small to constrain this 
scenario, the old GCs in NGC 5128 provide strong evidence 
for R— [3.6] colour bimodality, which can safely be inter- 
preted as metallicity bimodality. This confirms the spectro- 
scopic analysis of NGC 5128 GCs by Bcaslcy ct al. (200l). 
Strader et al. (2007) and Kundu & Zepf (2007) provide fur- 
ther evidence for GC system metallicity bimodality in mas- 
sive ellipticals, suggesting the Yoon et al. (2006) scenario 
does not apply to extragalactic (nor Galactic) GC systems. 

To conclude, despite a low IRAC pixel resolution and 
relatively short exposure times (e.g. 240s - NGC 4594, 72s - 
NGC 5128), the present work demonstrates that interesting 
science in the field of extragalactic GC system astronomy 
is indeed possible with Spitzer IRAC imaging. The largest 



prerequisite for such work is a good optical GC catalogue, 
because Spitzer photometry alone cannot be used for con- 
tamination removal for such faint objects. 
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